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16. Abstract 

In this study, a multi-layer emission model is described and 
applied to the measured emission from snow in an attempt to relate 
the absorption and scattering parameters of the snow medium to its 
density and wetness. 

The measured emission and ground-truth data used in this anal- 
ysis were those collected during a 24-hour diurnal experiment conduc- 
ted by Stiles and Ulaby (1980) on 2/17-2/18/77 in Steamboat Springs, 
Colorado. The emis.ion measurements were acquired at 37 GHz and 
H-polarization using a microwave radiometer mounted atop a truck- 
mounted boom. 

The Tinqa et al. (1973) mixing formula is used to calculate the 
dielectric constant of snow at 37 GHz. These values of the dielec- 
tric constant of snow along with the multi-layer emission model are 
used to estimate the absorption and scattering coefficients of snow 
at 37 GHz and their dependence on wetness. It is found that the scat- 
tering coefficient is comparable in value to the absorption coeffi- 
cient for dry snow. However, the absorption coefficient increases 
linearly with increasing snow wetness while the scattering (contd.) 
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Abstract (contd.) 


coefficient decreases linearly with increasing wetness. 

The estimated values of the scattering and absorption coef- 
ficients are then used to calculate the emission from each layer 
of the snowpack throughout the diurnal experiment. It is shown 
that for dry snow, the ground underneath the snowpack contributes 
about 45% of all measured emission while the rest is due to emis- 
sion from all the layers within the snowpack. However, when the 
wetness of the top 5 cm layer of the snowpack increases to 2% by 
volume, this top 5 cm snowlayer contributes more than 90% of all 
the measured eraissior. 
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^ MODEL DESCRIBING THE MICROWAVE EMISSION 
FROM A MULTI-LAYER SNOWPACK AT 37 GHz 


M. Abdel r«z1k 
Remote Sensing Laboratory 
University of Kansas Center for Research, Inc. 
Lawrence, Kansas 66045 


ABSTRACT 


In this study a multi-layer emission model is described and applied 
to the measured emission from snow in an attempt to relate the absorption 
and scattering parameters of the snow medium to its density and wetness. 

The measured emission and ground truth data used in this analysis 
were those collected during a 24-hour diurnal experiment conducted by 
Stiles and Ulaby (1980) on 2/17 - 2/18/77 in Steamboat Springs, Colorado. 
The emission measurements were acquired at 37 GHz and H-polarization 
using a microwave radiometer mounted atop a truck -mounted boom. 

The Tinga et al. (1973) mixing formula is used to calculate the 
dielectric constant of snow at 37 GHz. These values of the dielectric 
constant of snow along with the multi-layer emission model are used to 
estimate the absorption and scattering coefficients of snow at 37 GHz 
and their dependence on wetness. It is found that the scattering coeffi- 
cient is comparable in value to the absorption coefficient for dry snow. 
However, the absorption coefficient increases linearly with increasing 
snow wetneus while the scattering coefficient decreases linearly with 
increasing wetness. 



The estimated values of the scattering and absorption coefficients 
are then used to calculate the emission from each layer of the snow- 
pack throughout the diurnal experiment. It is shown that for dry snow, 
the ground underneath the snowpack contributes about 45% of all measured 
emission while the rest is due to emission from all the layers within 
the snowpack. However, when the wetness of the top 5 cm layer of the 
snowpack increases to 2% by volume, this top 5 cm snowlayer contributes 
more than 90% of all the measured emission. 
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1 .0 INTRODUCTION 


During the Winter of 1977, an experimental program was conducted 
at a test site near Steamboat Springs, Colorado Involving microwave 
measurements of snowpacks (Stiles and Ulaby, 1980). Among the Inves- 
tigat1on9conducted was an experiment consisting of a 37 GHz passive 
microwave measurement over a 24 hour diurnal period. In conjunction 
with the microwave measurements, several snow parameters were recorded 


which Include snow wetness, snow thermometric temperature and snow 
density at several depths in the snowpack. 

The purpose of this study Is to apply a multi-layer emission 
model to the measured emission from snow In an atteim)t to relate the 
absorption and scattering parameters of the snow medium to density 


and wetness. 


Chapter ^ covers a brief description of the diurnal experiment 
reported by Stiles and Ulaby (1980) and reviews relevant results 


obtained In other parts of the same Investigation. 

The development and applications of the microwave emission 
model are given In Chapter 4. Results estimated by the different 


models for the scattering and absorption coefficients of snow at 
37 GHz are shown. 


In Chapter 5, an evaluation of the different models which 
were developed is shown. A comparison between scattering and ab- 
sorption coefficients estimated in this study and those reported 
in the literature is also provided. 

Finally Chapter 6 summarizes the findings of this study and 
the errors involved in the development and application of the model. 
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2.0 DIELECTRIC PROPERTIES OF SNOW 


Snow is a mixture of ice» water and air. The absorption, emission 
and scattering of an electromagnetic wave by snow are governed by the 
geometrical and electrical (dielectric) properties of the snow medium. 
The relative dielectric constant of snow, K^, is in general a function 
of: 

(a) Microwave frequency; 

(b) Percentage volume of ice; 

'c) Percentage volume of air; 

(d) Percentage volume of free water; 

(e) Snow thermometric temperature; 

(f) Crystal size and structure of the snow medium; and 

(g) Presence of imourities. 

As will be seen in the next sections, the percentage of free water 
in a volume of snow has the strongest effect of all of the above para- 
meters on the value of K^. The dielectric properties of water and 
ice which constitute, along with air, the snow medium will be examined 
first, followed by an examination of the dielectric and attenuation 
properties of dry and wet snow at 37 GHz. 

2 . 1 Dielectric Properties of Water 

The dieiectric constant of water in the microwave region is 
governed by the Debye equation; 


- K 0 

s,' It. * j— 

1 - j2lTTf 2 itKq f 


( 2 - 1 ) 


1 



where 


= Complex relative dielectric constant of water, = 

* Optical limit of dielectric constant 
*^dc * Static limit of dielectric constant 
T = Relaxation time of the water molecule dipole 
= Permittivity of free soace 

0 

0 = Ionic conductivity 
f = frequency in Hertz 

The real and imaginary parts resulting from the previous equation are 


KJ, = K, + — ^ 

1 + (2-xf)^ 

(2nrf)(K. -KJ o 

K" = — 5 + 

w 1 . /o _r\^ 


1 + (2xxf)‘ 


2xKo f 


(2-2a) 

(2-2b) 


The Debye equation shows a single resonance phenomena at the relax- 
ation frequency fo given by: 


fo 


1 

7 ^ 


(2-3) 


The relaxation frequency is temperature dependent as shown in Table 2-1 
and is located in the microwave region. The limiting values on the 
equations for dielectric constant are 
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TABLE 2-1 

Relaxation Frequency of Water fo at 
Different Temperatures (Royer, 1973) 


T(°C) 1 

: folOHz) 

0 

8.5 

10 

11.7 

20 

15.8 

30 

21.2 

40 

27.0 

50 

! 

33.9 







for f«fo (2-4a) 

K " ^'^dc * 

= K^ for f»fo (2-5a) 

K" = (Kj^ - K^)/(2iTTf) for f»fo (2-5b) 


Figures 2.1 and 2.2 show the real and imaginary parts of the dielectric 

constant of water as a function of frequency and for different water 

tenroeratures. These curves were calculated by Royer (1973) using 

regression fits to obtain the values of K' and K" for several experi- 

mer.ters' data. The rate of attenuation in water is shown in Figure 2-3. 

It is seen that water is a very lossy m^^ium throughout the microwave 

region. This loss is due to the fact that the relaxation frequency 

(fo) is located in the microwave region. The attenuation of water 

4 4 

at 37 GHz is between 1.5 x 10 dB/meter and 2.0 x 10 dB/meter depend- 
ing on the temperature of the water. Because the dielectric constant 
of water is much larger than those of ice and air, small amounts of 
free water appearing in a snow medium can drastically alter the real 
and imaginary parts of the dielectric constant of snow and therefore 
change its emission, scatter and attenuation behavior. 

2.2 Dielectric Properties of Ice 

The dielectric constant of ice, like water, is governed by the 
Debye Equation 2-1. However the water molecules in ice are bound, 
and therefore not as mobile as those of free water. This results 
in a relaxation frequency for ice that is much lower than that of water 
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Relative Permittivity of Water at T = 0° C and 
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Figjre 2-? Re.ative Pfirmittivity of Water at T = 10° C and 
T = 30° C Using Debye Equation (Royer, 1973). 
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and is of the order of 10 KHz (Evans, 1965). The optical limit is 
therefore easily satisfied for the Debye equation and the real and imaginary 
parts of the dielectric constant of ice are given by equations (2-5a) 
and (2-5b). 

Measurements up to 183 GHz indicate that the real part of the 
dielectric constant of ice is approximately 3.17 in the microwave 
region and is independent of frequency and temperature as shown in 
Figure 2-4 (Evans, 1965). The imaginary part was found to bp tanpera- 
ture dependent at a frequency of 10 GHz showing a sharp in e near 
0®C (Lamb and Turney, 1949), Measurements of the dielectric ^unstant 
of ice are shown in Figure 2-4. Above 24 GHz, the dielectric constant 
measurements for ice are limited to the results of Perry and Straiton 
(1972). Table 2-2 shows their results at 35.3 GHz, These data are sus- 
pected to be in error (Gough, 1972) since they are anomalously different 
from both lower and higher frequency data. However Hallikainea (1977) 
used the results shown for the loss tangent in Figure 2-4 to calculate 
an average value for the loss factor of ice as a function of frequency. 

His results are shown in Figure 2-5. 

The rate of attenuation of tap-water ice in dB for different 
ice thickness at 35.3 GHz was also reported by Perry and Straiton 
(1972) and is shown in Figure 2-6. 

2.3 Dielectric Mixing Formulas 

Since snow is a mixture of ice, water and air, the varying 
volumetric properties of these components in snow make the dielectric 
properties of snow difficult to examine. In addition, the structure 
of snow is a significant factor in determining its dielectric properties. 
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Loss tangent of Ice versus radio fr^eocy. The quantity plotted 
vertically Is «) ^ frequency In iic/sec. 

On the hlqh- frequency tall of a relaxation spectrum this quantity 
Is constant: It has the further useful property that the 

attenuation of a radio wave (measured In dB/m) passing thiough 
the medium Is directly proportional to f tan 6. Temperatures a*'e 
marked in *C. 

L; Lamb (1946) and Lamb and Turney (1949) -5*C at low frequencies, 

0® to -190*C at high frequencies: distilled water. 

C: CuRining (1952) Distilled water, tap water, and melted snow 

(no observable difference) 

A: Auty ^nd Cde (1952) Conductivity water. Ice free fr^ stress. 
Llml #9 values plotted arbitrarily at 1,000 tim** the relaxation 
frequency. 

V: Von Hippel (1954) Conductivity w^ter. Ice not available, 

Y: Yashino (1961) Antarctic Ice core samples, not annealed, 
density 0.91g/cm3. 

U: Uestphal (private comoiunlratlon) Greenland Ice, annealed, 
density 0.90g/cm". 

B: Blue (1979) 


Figure 2-4 Dielectric Properties of Ice (Evans, 1965). 
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TABLE 2-2 

Relative Dielectric Constants of Ice at 35.3 GHz 
as Reported by Retry and Straiton, (1972) 



1.91 ±0.03 <4x10 

1.89 ±0.03 1.14x10 


1 

! tan 6 

Ice Source 

<2.1 X 10'^ 

6.03 X 10‘^ 

Deionized H 2 O 

H,0 
ap 2 







Figure 2-5 



Loss Factor of Fresh Water Ice as a Function of 
Frequency (Hallikainen, 1977). 
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Figure 2-6 Attenuation Curve (Tap-Water Ice) at 
35.3 GHz (Perry and Straiton, 1972). 
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Table 2-3 shows several mixing formulas that have been proposed to 
calculate the dielectric constant of a mixture consisting of two com- 
ponents (Poe, 1971; Sweeny and Colbeck, 1974; Tinga et al., 1973). 

The form number F, used in the first five formulas accounts for 
structure in the dielectric medium. The numerical values of the form 

number can vary from F = 0 (vertical particles), to F * 2 (spherical 
i*. 

partic^)'s), to infinity (for elongated horizontal particles). Values 
of F for freshly fallen snow are normally between 2 and 10 (Evans, 

1965). Edgerton et al . (1971) used a form number of F = 3? to get 
the best fit for their data for wet snow using the Weiner mixing formula. 

The mixing formulas generally show no dependence on frequency^ 
however, as the wavelength of interest approaches the order of the 
snow crystal size, scattering effects may become significant and the 
form number may require different values as a function of frequency. 

The wavelength at the desired frequency of 37 GHz is 0.81 cm which 
is of the order of snow crystal size; scattering, therefore, may 
become significant. 

The first five mixing formulas were used to fit dielectric 
constant data for snow below 10 GHz (Evans, 1965) and were reported 
to have worked satisfactorily. However, since few measurements of 
snow dielectric constant were conducted above 10 GHz, the mixing formulas 
were never used to fit any data for wet or dry snow at or near the desired 
frequency of 37 GHz. 

The Tinga et al. mixing formula was used to determine the dielec- 
tric constant of snow from the dielectric properties of ice and water 
at 5 GHz and 37 GHz (Tiuri and Schultz, 1980). The calculated values 
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TABLE 2-3 


Mixing Formulas 


Weiner (!'oe, 1971) 

Kc - 1 k, - 1 - 1 

k .. + F ■ ^i K. + F ” ■ k + F 

OS 1 3 


where 

p. = Volumetric Fraction of Ice 

1 

k^ = Dielectric constant of air 
k^. = Dielectric constant of ice 
*^ds ~ Dielectric constant of dry snow 
F = Form number 


Bottcher (Poe, 1971) 


3 k 


'‘i - '‘o 


ds ■ “i Ki * 2 kds 


Wet Snow (Weiner) (Poe, 1971) 


ws 


"ds 
'u * 0" 


(1 - 


w 




where 


U = 



* F 
TT 


k „5 = Dielectric constant of wet snow 

k = Dielectric constant of bulk water 
w 

p., = Volumetric fraction of water 

W 
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Pierce (Poe, 1971) 


^ds 


- k. + 


(1 - P^)(l 


F) 


0 - Pi)F 


(K - •«•) 
0 1 ' 


deLoor (Sweeny and Col beck, 1974) 



where 


k = Dielectric constant of bound water 
n 

A. = Particle Depolarization factor for each 
^ of the three major axes (Sweeney and Col beck, 
1974) 


Tinga et al . (1973) 

r/R 

- 2(y (V'HVS' -ft) 
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where 


k = Dielectric constant of water 
w 

k^ = Dielectric constant of ice 

R. = Radius of ice particle 

R = Radius of water covered ice particle 

= Weight percentage of liquid water contained 
in snow 

Pg = Density of snow 
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of the dielectric constant of snow were used to calculate the brightness 
temperature of snow using a radiative transfer model and assuming the 
snow to consist of spherical ice particles in air for dry snow. Wet 
snow was assumed to consist of the same ice spheres with a thin water 
shell surrounding them. Figure 2-7 shows the calculated dielectric 
constants of snow for the different snow wetnesses. The density and 
temperature of the snow were measured but the wetness and particle 
size have been selected for best calculated fit. It can be seen that 
a relatively good fit has been achieved for all cases. 

2.4 Dielectric Constant of Dry Snow 

Dry snow is a mixture of ice and air. Therefore, its response 
to temperature and frequency is similar to that of ice. However if 
the wavelength is of the order of the snow crystal size, scattering 
may cause a frequency dependence in snow. Due to the difficulty associ- 
ated with measuring the dielectric constant of snow while maintaining 
the environment constant during the measurement, only one set of dielec- 
tric constant measurements of snow above 10 GHz is reported. Table 
2-4 shows the results of the experiments conducted by Edgerton et 
al. (1971) at 37 GHz for dry snow. These measurements were made using 
an ellipsometer to measure the reflection coefficient. 

2.5 Dielectric Constant of Wet Snow 

Due to the !iiqh values of the dielectric constant of water, small 
amounts of water can significantly change the dielectric constant of 
snow. The lack of a good method to measure the wetness of snow along 
with the difficulty in keeoinq the environment constant during measurements 
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angle from nadir 



Thickness 

Density 

1fi«tness 

Ter^perature 

Particle 

radius 

Optical 

depth 

Diel.Re. 

Oiel.Ifli. 

Albedo 

(o) 

(k*/n>) 

(t) 

(C) 

(m) 

4.8/36.8 

4.S/36.8 

4.8/36.8 

4.S/36.8 

0.2 

350 

0.28 

-O.l 

O.S 

0.17/6.30 

1.S9/1.S7 

0.011/0/J064 

0.010/0.37 

0.18 

410 

1.1 

“0.5 

0.6 

0.66/13.1 

1.78/1. TO 

0.048/0.029 

0.0050/0.76 


Thickness 

M 

Density 

(Kg/m*) 

hetness 

(t) 

Teerperature 

(C) 

Particle 

radius 

(wa) 

Optical 

depth 

4.8/36.8 

Diel.Re. 

4.6/36.8 

Oiel . Im. 
4.8/36.8 

AltMdo 

4.8/36.8 

O.Ol 

350 

10 

0 

O.S 

0.10/1. SO 

2.07/1.67 

0.15/0.19 

0. 0027/0. T6 

0.19 

350 

6 

0 

0.5 

1.74/19.9 

1.92/1.62 

0.13/0.12 

0.0023/0.31 

0.18 

410 

2 

-0.3 

0.6 

1.07/15.8 

1.86/1.71 

O.03O/0.0S2 

0.0044/0.6S 


Figure 2-7 Measured and Calculated Brightness Temperatures of a 

?8 cm Thick Snow Field in the Morning and in the After- 
noon. The Density and Temperature of Snow Were Known. 
The Wetness anu Particle Size have been Selected for 
Best Fit (Tiuri and Schultz, 1980). 
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TABLE 2-4 


Dielectric Constant of Dry Snow at 
37 GHz as Reported by 
Edgerton et al . (1971 ) 



k" 

s 

Temperature 
° C 

niii^^^iiiiin 

Grain Size (mm) 

2.77* 

< 0.03 

-2 

0.5 

1 ,2 

1.9 

< 0.03 

1 

-10 

0.5 

0.5 

. . i 


•k 

The high values were explained as resulting from difficulty in sample 
preparation. 
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make dielectric constant measurements of wet snow difficult obtain. 
Measurements of dielectric constant of wet snow above 10 GHz are nonexis- 
tent. Therefore, there is no data for the dielectric constant of snow 
of any wetness at or near the desired frequency of 37 GHz. 

2.6 Attenuation Through Snow 

The absorption coefficient for a lossy medium is given by: 



However, experimental evidence indicates that for wet and dry snow: 


k;«k 


f 

S 


This reduces equation (2-6) to: 



where 

a. * Field absorption coefficient, nep/m 

O 

X * Wavelength (meters) 

= k; - j k; 

= Complex dielectric constant of snow 
Attenuation rates in dB/m are given by: 



(2-7) 


( 2 - 8 ) 


(2-9) 
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Power loss measurements made through a layer of snow are composed of 
two parts, mismatch loss and attenuation loss. If the power loss 
measurentents are conducted with coherent transmission, the two types 
of loss cannot be seoarated because of the interference effects caused 
by multi ole reflections. 

Linlor (1980) measure"! the rate of attenuation in dB/cm as a 
function of snow wetness and at selected frequencies between 4 and 
12 GHz using phase shift and transmission loss measurements through 
two identical containers. One of the containers was empty and used 
as the reference unit and the other was filled with the snow sample. 
Figures 2-8 and 2-9 show the results of his experiment. These atten- 
uation curves were plotted using an empirical relation derived by 
Linlor (1980) ano given as: 


where 


Attenuation(d3/cm) = 


0.045(f-4) + 0.066 


= Percentage wetness by volume 
f = frequency in GHz 


(2-10) 


Notice the increase in attenuation with increasing snow wetness and 
increasing frequency. Figure 2-8 shows that the attenuation through 
dry snow (zero wetness) is zero. This, according to the author, is 
not the case but the attenuation was less than 1 d3 and so could not 
be measured using the available measuring instruments. 

The Tinga et al. (1973) mixing formula along with Equation (2-8) 
are used to calculate the attenuation through snow as a function of its 

3 

wetness at 37 GHz, 18 GHz and 12 GHz, for snow with a density of 0.21 q/cm , 
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Figure 2-8 


Figure 2-9 



SNOW A'ETNESS (volume %) 


Variation of Attenuation with Snow Wetness at Selected 
Frequencies (Linlor, 1980). 



Variation of Attenuation with Frequency at Selected 
Snow Wetness (Linlor, 1980). 
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Figure 2-10 shows the results of this calculation. As seen by the figure, 
there is an increase in attenuation with increasing wetness of the 
snow and with increasing frequency. However, at 12 GHz, the attenua- 
tion starts to saturate at a percentage wetness, M^, of 3%. In 
Figure 2-11, the attenuation in Nep/cm is plotted as a function of 
at 12 GHz using both Equation (2-10) of Linlor (1980) and the Tinga 
et al. (1973) mixing formula. It is seen that the attenuation values 
obtained by the mixing formula are higher in value than those obtained 
by Linlor's empirical formula for percentage wetness by volume of 
snow that are less than = 2.75. 

Battles and Crane (1965, 1965) used interferometry techniques 
to measure loss from artificially created snow. Figure 2-12 shows 
their results in dB/ft as a function of snow temperature at 35.2 
GHz. The loss shows a slow increase with increasing temperature 
until the temperature nears the melting point. Near 32° F, the 
melting point of ice, a sharp increase in loss through snow is seen 
which is due to the appearance of moisture in the sr.ow. Table 2-5 shows 
the reseats of another experiment reported by Battles and Crane (1966), 
in which loss from different types of snow is reported at 35.26 GHz. 

The higher losses for locally packed and ice crystal cases were explained 
to be the result of scattering of the wave by the snow crystals. 

Currie et al . (1977) made measurements using a pulsed radar 
operating at 35 GHz. Their loss measurements were calculated by 
comparing the return from a corner reflector placed below the snow 
layer with the return measured from the same corner reflector placed 
above the snow layer. Table 2-6 shows their results. The loss measure- 
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Figure 2-1 
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Mv(« 


1 The Attenuation Through Snow as a Function of its Wetness 
as Calculated by Linlor's (1980) Empirical Formula and 
by the Tinga et al, (1973) Mixing Formula at 12 GHz 
and Snow Density of 0.21 g/cm^. 
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-12 Absorption of Radiation by Snow as a Function of Temper- 
ature (Battles and Crane, 1965). 
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TABLE 2-5 


Loss Caused by Snow and Ice at 35.26 GHz 
(Battles and Crane, 1966) 


Type 

Thickness 

(cm) 

Density 

g/cm^ 

Loss 

(dB) 

Fine, loose snow 

14.0 

0.2 

2.5 

Locally packed loose snow 

14.0 

0.33 

7.2 

Large ice crystals 

14.0 

0.39 

9.2 

Packed snow 

14.0 

0.47 

1.4 

Ice 

5.1 

0.92 

1.2 
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TABLE 2-6 


Loss Measurements at 35 GHz 
(Currie et al., 1977) 


Snow Condition 

Layer Thickness (cm) 

Loss (dB) 

Wet Snow Crust 

3.5 

11 

Dry Snow Crust 

4.6 

21 

Dry Snow Crust 

5.0 

5 

Dry Snow Crust 

i 

3.5 

3 
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ments are due to both mismatch and attenuation losses and, therefore, 
an attenuation coefficient cannot be calculated. 

Stiles and Ulaby (1980) measured path loss at 35 GHz as a func- 
tion of depth for three snow conditions using transmission techniques. 
Figure 2-13 shows their results In dB as a function of snow depth. 

Stiles and Ulaby (1980) also calculated the emissivlty of snow 
at 37 GHz as a function of the angle of Incidence of the radiometer 
relative to nadir. The emissivlty e(e) was calculated using the equation 
below 


e(9) = 




( 2 - 1 1 ) 


where 

Tap = Apparent temperature of snow measured by radiometer; 
°K 

T^l^y = Sky brightness temperature; °K 
Tq = Physical temperature of snow; °K 

Stiles and Ulaby (1980) developed a simple emission model for a homo- 
genous layer of snow and t!ie ground underneath It, and were able to 
fit their emissivlty data calculated using the previous equation to 
the following equation developed by their model for emissivlty as 
a function of the angle of incidence e: 
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Total Loss (cIB) 





where A and B represent the constant terms inside the parenthesis and 
T,, * Snow-air power transmission coefficient 

Sc* 

Tyj = Ground-snow power transmission coefficient 

K - Snow absorption coefficient (nep/cm) 

a 

Kg = Snow extinction coefficient (nep/cm) 

' ("a * S> 

-1 2 

Kg = Mass extinction coefficient of snow (nep g cm ) 

K^ = Snow scattering coefficient (nep/cm) 

W = Snow water equivalent, cm 

0 = Angle of incidence, relative to nadir (degrees). 

Figure 2-14 shows a plot of measured radiometric emissivity for dry 
snow at 37 GHz and for e = 27° and 57° as a function of the water 
equivalent W of snow. The best fit curves obtained using Equation 
(2-12) to fit the emissivity data are also shown. The values of A 
and B estimated by regression analysis to produce the best fits are 
shown in Figure 2-14 and Table 2-7. From these values it is seen 
that at 9 = 27°: 


^ = 0.517 (2-13a) 

sa Kg 


- T,, — = 0.481 
sa gs sa <g 


(2-13b) 


sec 27° = 0.0235 (2-13c) 

Assuming t = 0.99, for the snow density measured during the experi- 

3 

ment (0.41 g/cm ), the above three equations are used to calculate k , 

a 
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Emissivity 



Water Equivalent W (cm) 

Figure 2-14 Measured Radiometric Emissivity Response to Dry Snow Water 
Equivalent at 37 GHz (Stiles and Ulaby, 1980). 
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TABLE 2-7 


Values of A, B, ice» '<^a •'s Generated by Fitting 

the Experimental Values of e at 37 GHz to Equation 2-11; 
(Stiles and Ulaby, 1980) 



■'e 

nep/cm 


Maximum 



27° 0.517 0.481 8.60 x 10’^ 4.50 x 10’^ 4.10 x 10”^ 0.02 

57° 0.586 0.273 1.38 x 10’^ 8.17 x 10“^ 5.63 x 10"^ 0.02 


*Maximum Error = |e^ - eni||; = Calculated from best fit eq. 2-12. 

= Measured (equation 2-11). 
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< and K for snow at 9 = 27°. The results of the calculations are 
s e 

shown in Table 2-7. Similar calculations are carried out for 9 = 57° 
and the values obtained are also shown in Table 2-7. 

Attenuation rates for snow were calculated using the temperature 
measurements of a snow layer of known thickness over a metal plate 
(Hofer and Matzler, 1980). The attenuation rates were calculated 
using two simple models proposed by Hofer and Matzler (1980). The first 
model included the effects of multiple reflections but neglected inter- 
ference effects, whereas the second model was a radiative transfer 
model which included volume scattering -c and absorption < Figure 
2-15 and 2-16 show the results. <2 is the power damping coefficient 
which combines the effects of absorption and scattering and is given 
by (Hofer and Matzler, 1980): 



is the power attenuation coefficient calculated using the first 
model. Hofer and Matzler (1980) reported moist snow with wetness 
of 1 to 3 percent by volume exhibited power damping coefficients 
greater than 30 dB/meter at 36 GHz. 
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4.9 

Figure 2-1 


.5 21 36 

Scattering, Absorption and Damping Coefficients for 
Dry Winter Snow (Hofer and Matzler, 1980). 
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3.0 EXPERIMENT DESCRIPTION AND RESULTS 


This chapter provides a description of the microwave diurnal 
experiment conducted by Stiles and Ulaby (1980) on 2/17 - 2/18/77 
and on which this study is based. The first section gives a brief 
description of the sensor used, test site and the different microwave 
and ground truth experiments conducted. Later sections of the chapter 
illustrate the results of the experiment. 

3 . 1 Description of the Microwave Radiometer and Ground Truth Data 

A microwave radiometer was mounted atop a truck-mounted boom and 
ured to acquire ^ppp/lent temperature (T ) data at a test site near 
Steamboat Springs, Colorado. The radiometer used during the experi- 
ment was a dual -polarized Dicke type manufactured by Aerojet General 
Corporation. It operated at 37 GHz, had a bandwidth of 300 MHz, 
a sensitivity of 0.5°K, an accuracy of 1°K and operated over a temper- 
ature range of 0-500°K. The radiometer had an approximate gain of 
0.010 volts/K and used an automatic gain corarol mechanism. Calibra- 
tion of the radiometer was conducted at the University of Kansas after 
completion of the experiment (Stiles and Ulaby, 30). 

The test site chosen for the experiment was a 40-acre hayfield 
near Steamboat Springs. The surface of the field was very flat and 
was covered with close-cut hay, approximately 6 cm in height. A snow- 
pack approximately 30 cm in depth lay on the surface of the field 
during the 2/17 - 2/18/77 diurnal experiment. 

The diurnal experiment consisted of continuous data acquisition 
over a 28-hour period commencing at 6:00 AM on 2/17//7. T data 

u yj 
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was acquired at 37 GKz at both the horizontal and vertical po1ariza> 
tions and at 0°, 20® and 50® angles of incidence. Ground truth 
data were gathered in a snow pit on the test site and included measure- 
ments of snwn density profiles, snow wetness profiles, snow temperature 
profiles along with soil moisture samples and soil tenqierature in 
“■y jO 5 cm of the soil. These measurements were obtained at 
hourly intervals except for the snow wetness measurements which required 
a slightly longer duration. Snow density profiles were obtained by 
taking horizontal cores from each layer in the snowpack. Snow wet- 
ness was measured wiw.i a freezing calorimeter at 5 cm intervals ver- 
tically throughout the snowpack. Snow ten^erature was measured 
using a digital thermometer and was monitored at 2 cm intervals ver- 
tically throughout the snowpack. Since water equivalent is the total 
amount of water contained in the form of snow per unit area, it is 
calculated as the product of density and depth of each snow layer. 
Results of T^p and ground truth measurements are discussed below 
(Stiles and Ulaby, 1980). 

3.2 Diurnal Variation of at 37 GHz 
afi 

Since the model developed later in this paper uses only the T^p 

data obtained at 37 GHz in the he izontal polarization mode, only t’le 

horizontal polarization data are illuitrated here. Figure 3-1 shows 

the apparent temperature T^p response at 9 = 0®, and 50® throughout 

the diurnal experiment. Notice the increase in T starting at 

ap 

1000 hours and ending at 2200 hours. This increase in T,^ is due to 

ap 

the appearance of water in the snowpack as a result of warmer snow 
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c-imjre 3“1 Apparent Temperature as a Function of Time at 37 GHz, H Polarization and for 
« no. 2(P and 50O (Stiles and Ulaby, 1980). 





temperatures. Notice also the little difference between the T^p response 
at e - 0^ and 20^. Figure 3.2 shows the snow wetness in individual 
5 cm layers in the top 15 cm of snow as a function of time. By com- 
paring Figures 3-1 and 3-2 one notices that T^p is shown to lag behind 
N^(0-5 cm), the percentage wetness by volume of the snow in the top 
5 cm of the snowpack, by about two hours. T^p exhibits sharp increases 
and decreases as the snow starts to melt and refreeze respectively. 

3.3 Ground Truth Diurnal Variations 

As mentioned in the section above. Figure 3-2 shows the wetness 
of the snowpack at the depths of 0-5 cm, 5-10 cm and 10-15 cm 
from the snow surface. The snowpack at depths lower than 15 cm was 
frozen throughout the experiment and showed no wetness. is 
defined as the volume percentage of liquid water contained in a 
unit volume of snow. 

Snow thermometric temperatures measured every 2 cm from the top 
of the snowpack to a depth of 5 cm below the soil surface were used 
to generate the average values shown in Table 3-1. In this table, 
the snow therraometric temperatures are averaged to give an average 
temperature for the snow every 5 cm from the snow surface. Table 3-1 
shows 21 different data sets representing snow thermometric temperatures 
throughout the diurnal experiment. Figures 3-3, 3-4, 3-5 and 3-6 
show the snow thermometric temperature along with the snow wetness 
at intervals of 5 cm of snow and at the hours of 5:30, 10:00, 18:00 
and 6:45 respectively. Notice that as the snow starts getting warm 
and starts melting, it warms up first at the surface layer of snow 
and then the lower layers of the snow begin to warm up. Later in 
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Figure 3-2 Snow Wetness My as a Function of Time in the Top 15 cm of Snow (Stiles and Ulaby, 


TABLE 3-1 


Snow Therraometric Temperature in °K for Different Depths 
Within t\r ' Snowpack on 2/17/77 - 2/18/77. 

(-uiles «rd Ulaby. 1980) 


Time 

0-5 cm 

5-10 cm 

10-15 cm 

15-30 cm 

0530 

267.9 

267.2 

268.1 

269.3 


270.8 

269.5 

269.5 

270.8 

1000 

273.8 

272.5 

271.4 

272.0 

1100 

274.0 

272.3 

271.8 

272.3 

1200 

274.0 

273.3 

272.6 

272.6 

1300 

1 

274.0 

273.3 

272.8 

272.8 

1400 

274.1 

273.4 

273.2 

272.9 

1500 

274.0 

273.4 

273.3 

272.8 

1600 

274.0 

273.4 

273.3 

273.0 

1700 

273.7 

273.4 

273.2 

273.0 

1800 

273.4 

273.3 


273.0 

1900 

273.2 

273.2 

273.1 

273.0 

2000 

273.1 

273.1 

273.1 

273.0 

2100 

272.0 

273.0 

273.0 

273.0 

2200 

271.1 

273.0 

273.0 

272.9 

0100 

267.8 

270.1 

271.9 

272.6 


262.3 

264.4 

268.0 

271.2 

0615 

260.0 

263.6 

267.3 

270.5 

0645 

259.8 

263.5 

267.3 

270.4 

0745 

262.9 

262.7 

265.9 

269.9 

0820 

266.1 

1 264.4 

266.9 

269.7 
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Snow Temperature In ®C 



Figure 3-3 Snow Thermometric Temperature and Wetness at 5:30 AM 

on zmm. 
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Snow Temperature In ®C 



-5 30 25 20 15 10 5 0 


Depth of Snow in cm. 

Figure 3-4 Snow Thermoi.^tric Temperature and Wetness at 10:00 AM 
on 2/17/77. 
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Mv (%) 



Snow Temperature in ®C 



-5 30 25 20 15 10 5 0 


Depth of Snow in cm. 

Figure 3-5 Snow Thermometric Temperature and Wetness at 1800 Hours 
on 2/17/77. 
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3.0 

2.6 

2.2 

1.8 

1.4 

1.0 

0.6 

0.2 

0.0 


Figure 3-6 Snow Thermometric Temoerature and Wetness at 6:45 AM 
on 2/18/77. 
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the day as the surface layer becomes very wet, the lower layers warm 
up to almost the same temperature as the top 5 cm snow layer as shown in 
Figure 3-5. Similarly at the refreezing stage after the sun goes down 
and the snow starts drying up in moisture, the top 5 cm layer starts 
decreasing in temperature first followed by the lower snow layers. 
Therefore, the refreezing process starts, like the melting process, 
from the snow surface down to the lower layers of snow. 

Snow density was measured for each physical 
layer from the top of the snowpack down to the soil surface. The 

3 

average density was found to be 0.21 g/cm and did not change signif- 
icantly at different depths in the snowpack or with time throughout 
the diurnal experiment. Therefore, the water equivalent of the total 
snowpack did not change and was measured to be 6.3 cm throughout the 
entire experiirent. Also, the snow depth remained constant at 30 cm. 

3.4 Wetness Response 

Figures 3-7, 3-8 and 3-9 illustrate the response of T^^ to snow 
wetness of the top 5 cm layer at 9 = 0°, 20° and 50° respectively. T^p 
is shown to exhibit a hysteresis-like pattern. This hysteresis-like 
pattern is due to the fact that as the snow gets wet, the microwrve 
radiometer operating at 37 GHz becomes sensitive to a snow depth that 
is smaller than the top 5 cm layer of the snow from which the snow wetness 
measurement was obtained. The snow depth to which the sensor is sensi- 
tive when the snow is wet depends on the penetration depth of the wave 
in the wet snow layer. In other words, the snow wetness actually measured 
is the average wetness of the top 5 cm of snow. The wetness values for 
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Mv (0-5 cm) )?9 


Figure 3-7 Tgp Response to Snow Wetness M , Showing the Hyster- 
esis Effect at 0 = 0°. '' 
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Figure 3-8 Ta^ Response to Snow Wetness M . Showing the Hyster 
esTs Effect at = 20°. '' 
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Figure 3-9 Tap Response to Snow Wetness M , Showing the Hyster- 
esis Effect at 9 = 50°. ^ 
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the top 5 cm snow layer therefore gives a value of wetness that is 
stroller than the actual wetness of the snow depth ^ which the sensor 
is sensitive at the melting stage and a larger value of wetness et 
the refreezina stage. This causes the T._ response to lead in the 

ap V 

melting stage and to lag behind it in the refreezing stage which is 
shown in Figure 3-1 and is the cause of the hysteresis effect. 

3.5 Depth Response 

Stiles and Ulaby (1980) conducted an experiment in which they 

measured T^p at 37 GHz for different snow depths and different angles 

of incidence. The depth of snow was varied by piling dry snow to a 

depth of 170 cm. Apparent temperature T shows an exponential- like 

decrease with increasing water equivalenc, as shown in Figure 3-10. 

Data collected from an adjacent undisturbed snowpack shows values 

of T^^ that are considerably lower than those from the snow pile, 
ap 

This is due to the presence of ice layers in the natural snow and 
to differences in layerinn and in the distribution of crystal sizes 
which contribute to a larger scattering loss at 37 GHz and, therefore, 
lower emission. 
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Figure 3-10 Radiometric Apparent Ter ^rature Response to 
Sn; • Water Equivalent at o7 GHz (Stiles and 
1980). 
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4.0 MULTILAYER EMISSION MODEL DEVELOPMENT AND TESTING 


The relationshio between the measured microwave data (T ) 

ap 

from a snowpack and the physical and dielectric parameters of the 
snowpack can be examined through the development of models to explain 
this relationship. In the first section of this chapter, an emission 
model for snow is developed in a more general form. Later sections 
of this chapter show the use of a mixing formula to calculate the 
complex dielectric constant of snow at 37 GHz for different snow 
wetness. The calculated complex dielectric constant of snow is used 
along with the different forms of the emission model developed to 
estimate the absorption and scattering coefficients of snow at 37 GHz and 
their dependence on snow wetness. These coefficients are then used 
along with the snow model to estimate the values of T^^ measured through- 
out the diurnal experiment. 

4.1 Development of the Multilayer Emission Model for Snow 

Consider a snowpack of height d above the ground and consiscing 
of N different snow layers. It is assumed that: 

(a) Each snow layer in the snowpack has a uniform thickness d. 
cm, where i stand, for the snow layer number from the soil surface 

up to the snow surface. Dielectric and thermal properties of each 
snow layer are assumed to be constant within the layer, at a given ..ime. 

(b) The surface of the snow and the soil underneath the snovpack 
are smooth. 

(c) Multiple reflections at layer boundaries are ignored. 

(d) The radiation received by the radiometer is incoherent. 
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(e) Since the radiometer used In the experiment was a ground 
based system, there is neither attenuation nor emission between the 
surface of the snowpack and the radiometer. 

(f) The snowpack is a layered scattering volume with and 
Kj. taken to be the absorption and scattering coefficients of the i^^ 
snow layer respectively, tc^. is the sum of and and is the ex- 
tinction coefficient of the i^^ snow layer. Explicitly, diffuse scatter 
is ignored^ however it is understood that the above coefficients are ef- 
fective values that at least partially compensate for this approximation. 

The apparent temperature measured by the radiometer at different 
angles of incidence 9 is given bv: 

Tap(‘>) = Tb(o) * T^^O) (4-1) 

v*here Tg is the brightness temperature due to emission from the 
snowpack and the ground underneath it, and T^^ is the component of 
the downward emitted sky radiation scattered by the scene in the 
direction of the radiometer. The brightness temperature of the scene 
measured at an angle of incidence 0 is given by: 

Tb*”*) = ^^(N-) (4-2) 

where 

Tg(N^) = The brightness temperature measured directly 
above the snow-air boundary, 

= The power transmission coefficient from layer N, 
the uppermost layer of the snowpack to air 

Tg(N') = The brightness temperature measured directly below 
the snow-air boundary 
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Tg(N~) is the sum of the components due to emission by all the layers 
within the snowpack T^(SUM), and the comoonent due to the emission 
by the ground Tg, therefore 

Tg(N') = T^(SUM) + Tg (4>3) 

Self-emission by a thin layer of snow of thickness Ah within the 
i^^ layer of snow is given by: 

Tsel = T, j 1 - EXP P, sec 6| Pb )| 5 T, . p, Sec e: 

(4-4) 

The letter i above is an indication of the position of the snow layer 

considered within the snowpack, starting with i = 1 for the snow 

layer on the bottom of the snowpack directly above the ground surface, 

and i = N for the snow layer at the top of the snowpack directly below 

the snow-air boundary. T- is the physical temperature of the i^^ 

layer of snow, p. is its density, k^. is its absorption coefficient and 

el is the angle of propagation in the i^^ snow layer relative to nadir. 

The total emission by the i^^ snow layer thickness d^ is 

(Stiles and Ulaby, 1980): 

di 

° ’’i/’ f 'ai 
0 

(4-5) 

where is the scattering coefficient of the i^^ layer of snow. 

After integration. Equation (4-5) reduces to: 




dh 
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The radiation emitted by the i^^ snowlayer travels through the snow layers 
above it before it is seen directly below the air-snow boundary. The 
radiation emitted by the i^^ layer of snow is, therefore, attenuated 
by the factor 


Exp 


i * -sV i "j 'i 

j=i+1 


(4-7) 


The above factor is the loss due to the orooagatirn of the radiation 
emitted by the i^*' snow layer throuqh the (N-i) snow layers above it. 
Therefore 


1 N-1 


Ts(SUM) =V (Il'i(f.,)) ■ Exp 


i = l ' 

0 . d. Sec 


■ 

j 


S I ‘^’aj '"sj^ 

j = i+l 

(4-3) 


^i(i+l) power transmission coefficient across the boundary between 

the i^^ snow layer and the snow layer just above it. Similarly i^ is the 
ground thermoraetric temperature undern®atn the snowpacx and is the 
power transmission coefficient at the ground-snow interface, the ground 
contribution measured directly below the snow-air boundary is: 
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T. = 




N-1 

^9 


Exp 




Sec et 


1=1 


(4-9) 


Tg^ is the power transmission coefficient at the boundar; between the 

ground and the lowermost snow layer (layer =1). 

The value of T^. given by Equation (4-6) is substituted into 

t 

Equation (4-8), and this later equation givina an expression for T (S'JM) 

A S 

along with Equation (4-9) giving an expression for T^, are substituted 

into Equation (4-3) producing a complete expression for Tg(N'). The 

final expression for Tg(n") substituted into equation (4-2) to 


mve: 


1 /N-1 




1 - Exp 

i 






''j 'j| * ’na 'ql '4'!,''^ ('*'>) 


(4-10) 


S |-f-aiSl> -’i '*i ’i| 

‘i=l 

The above equation describes the emission from a snowpack ot N 
different layers above the ground and the emission from the ground underneath 
the snowpack. This equation when substituted into eauation (4-1) 
gives the general expression of the model developed. 

T^c> the component of the downv/ard emitted sky radiation scattered 
by the scene in the direction of the radiometer, shown in equation (4-1) 
is given as: 


?S 



(4-11) 


Tsc'‘> = ■ ’na> ^ky 

T^ky is the downward emitted sky radiation. Equation (4-11) when added 
to Equation (4-10) produces an expression for the apparent temperature. 
Tap, measured by the radiometer. 

4.2 Calculations of the Dielectric Constant and the Power Absor, on 
Coefficient of Snow 

In order to apply the ei.iission model developed in the previous 

section and described by Equation (4-1), calculations of the power 

transmission coefficients at every boundary connecting two different 

snow layers and at the snow-air interface and the ground-snow interface 

are neeued. In addition, values of the propagation angle in every 

snow layer and the absorption and scattering coefficients of every snow 

layer need to be calculated. To obtain all these values, the complex 

dielectric constant of snow at 37 GHz is needed as a function of snow 

wetness and snow density. Since measurements of the dielectric constant 

of snow at 37 GHz do not exist, one of the mixing formulas shown in 

Table 2-3 will be used. The Tinqa et al . (1973) mixing formula is 

used to obtain the dielectric constant of snow as a function of its 

wetness and density at 37 GHz. The Tinga et al., mixing formula 

is osec nere because it was shown to yield good results by Tiuri 

and Schultz (1980) when used at 37 GHz. The radius of the ice particle 

is taken to be 0.5 mm. the density of snow, was measured during 

the ground trutn experiment, and previously reported in Section 3-3. 

equal to 0.21 o/cm^. The value of k , the dielectric constant of water 

w 
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at 37 GHz, is obtained from Figure 2-1 and is k = 9.55 - j 19.10. 

w 

The value of k. , the dielectric constant of ice at 37 GHz, is obtained 
from Figure 2-5 an'* is equal to k. = 3.15 - j 0.003. The values of the 
dielectric constant of snow, obtained from the mixing formula, as a 
function of the percentage wetness by volume of the snow is shown in 
Table 4-1 and Figures 4.1 and 4.2. 

The imaginary part of the dielectric constant of snow, k“, shows 
a very small value when the snow is of zero wetness. However, this 
value increases substantially as the snow wetness increases by a very 
small amount, as shown in Table 4-1 for a percentage of snow wetness 
by volume of 0.1%. The difference in wetness between completely dry 
snow and that of 0.1% wetness by volume cannot be detected by any avail- 
able experimental methods. This leads us to suspect that the value of 
the loss factor calculated for completely dry snow is unrealistic and, 
therefore, should be replaced by a value comparable to that calculated 
for a snow wetness of 0.1% by volume. If a linear fit is used to fit 
the values of k^ at percentages of snow wetness by volume of 0.1% 
and 0.2% and then the fit is extended to oredict the value of k" 
for snow of zero wetness, the value of k^' predicted is 0.0015 as 
shown in Figure 4-3. This value is very close to that reported by 
Sweeny and Colbeck (1974) who experimentally measured the loss factor 
for dry snow at a frequency of 6 GHz and found it to be 0.003. Hence, 
the value k^ = 0.003 is the one chosen in the investigation for = 0. 

The values of k^ and k^ obtained by the mixing formula of Tinga 
et al. (1973) and modified at one value of k^ for dry snow are sub- 
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TABLE 4-1 


TKe Dielectric Constant of Snow 
as a Function of the Percentaqe Wetness by Volume of Snow 
at 37 GHz* Calculated by the Tinga et al. Kixing Formula, 


N 

V 



0.0 

1.317 

0.2844 X 10'^ 

0.1 

1.319 

0.5352 X 10*^ 

0.2 

1,321 

0.0113 

0.3 

1.323 

0.0168 

0,4 

1.326 

0.0221 

0.5 

1.328 

0.0273 

0.6 

1.331 

0.0325 

0.7 

1.334 

0.0376 

0.8 

1.337 

0.0426 

0.9 

1.339 

0.0475 

1.0 

1.342 

0.0523 

1.1 

1.344 

0,0547 

1.2 

1.345 

0.0571 

1.3 

1.348 

0.0617 

1.4 

1.351 

0.0662 

1.5 

1.355 

0.0707 

1.6 

1.358 

0.0751 

1.7 

1.361 

0.0793 

1.8 

1.364 

0.0835 

1.9 

1.368 

0.0876 

2.0 

1.371 

0.0916 

2.1 

1.375 

0.0954 

2.2 

00 

0.0992 

2.3 

1 .332 

0.1029 

2.4 

1.385 

0.1065 

2,5 

1.389 

0.1101 

2.6 

1.392 

0.1135 

2.7 

1.39u 

0.1168 

2.8 

1.400 

0.1200 

2.9 

1.403 

0.1232 
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Figure 4-2 


The Imaginary Part of the Dielectric Constant of Snow 
k^, as a Function of M^, the Percentage Wetness by Volume 
of Snow as Calculated by Tinga et al . (1973) Mixing 
Formula. 






stituted in equation (2-8) to give the values of a., the field absorption 
coefficient of snow as a function of the snow wetness. The corresoonding 
values of k., 

a 


'a ■ 

where k. is the power absorption coefficient in Nep/cm, are shown in 

a 

Table 4-2. Figure 4-4 shows a plot of k. in Nep/cm versus the 

a 

percentage wetness by volume of snow, M^. This figure shows that the 
dependence of k . on wetness is linear and its linear best fit is given by: 

a 

Kg = 0.0203 + 0.2741 , Nep/cm (4-13) 

The equation abcve is used throughout the rest of this investi- 
gation to describe the power absorption coefficient of snow at 37 GHz. 

4 . 3 Emission Model for a Single Homogeneous Snow Layer : 

During the diurnal experiment conducted on 2/17 - 2/18/77 by 

Stiles and Ulaby (1980) the total depth of the snowpack above the 

ground was 30 cm. The snowpack w<is divided into four layers. The 

thermometric temperatures of each snow layer were measured throughout 

the experiment and are shown in Table 3-1. The percentage wetness by 

volume, M^, of each snow layer .werer measured throughout the experiment 

and shown in Figure 3-2. T measurements for o = 0°, 20° and 

ap 

50°, as measured by the radiometer are reported in Figure 3-1. Com- 
bining the results shown in Table 3-1 and Figure 3-2, produces 21 complete 
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TABLE 4-2 


Values of in Nep/cm as a Function of Snow Wetness 

O 


{%) 

K, in Nep/cm 

a 

0.0 

0.2030 X 10'^ 

0.11 

0.3949 X 10"^ 

0.21 

0.7649 X 10'^ 

0.32 

0.1129 

0.42 

0.1487 

0.53 

0.1839 

0.63 

0.218‘ 

0.74 

0.2523 

0.84 

0.2855 

0.95 

0.3181 

1.05 

0.3500 

1.15 

0.3812 

1.26 

0.4117 

1.36 

0.4416 

1.47 

0.4707 

1.58 

0.4992 

1.68 

0.5269 

1.78 

0.5539 

1.89 

0.5803 

1.99 

0.6059 

j 2.10 

0.6309 

i 2.20 

1 

0.6551 

1 2.31 

0.6787 

2.41 

0.7016 

2.52 

0.7237 

1 

2.62 

1 ... ’ 

0.7453 j 

■ — — - i 
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iCa Nep/cm 



ground truth and microwave measurements at different times in the 
diurnal experiment. For eight out of the 21 different measurements, 
the snow was totally dry in all its layers within the snowpack. 

In this first application of the model described by equations 
(4-10) and (4-11), the snowpack is assimied to be made up of one honw- 
g 90 US snow layer that is 30 cm thick. The thermometric temperature 
of this homogeneous snow layer, its density and its wetness are assumed 
similar to those measured for the top 5 cm snow layer of the snowpack. 
In other words, an attempt will be made in this section to see whether 
the temperature, wetness and density of the top 5 cm layer of snow 
can be used to represent the whole snowpack. 

When Equation (4-10) is used to describe the emission from one 
layer of snow above the ground it reduces to: 


° sa '^1 ( j| * ^s|)|' I * ''s'* 

'“l “ijh "9S 'sa ^’'“[-(■'al * "sl> 

0 ^ d-| Sec e-j 1 


(4-14) 


where 

T_, = Power transmission coefficient at the snow-air bour.oary, 

Sa 

Tgs " Power transmission coefficient at the ground-snow boundary, 
' The snow thermometric temperature in °K. 


The equation above is solved in the following way: 
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(a) T^, the snow thermometric ten^erature values were shown In 
Table 3-1. values used are ttose reported for the top 5 cm snowlayer. 
These temperatures are assumed indicative of the whole snowpack. 

(b) pp the density of snow was measured during the ground truth 

3 

experiment and reported equal to 0.21 g/cm . This value did not change 
throughout tl . diurnal experiment. 

(c) dp che total depth of the snowpack was measured throughout 
the experiment and was equal to 30 cm. The depth of the snowpack did 
not change throughout the diurnal ex. erlment. 

(d) Tg, the therraometrlc temperature of the ground underneath 
the snow was measured for the top few centimelars of the ground and 
was equal to 273°K. That Is the ground was frozen throughout the 
diurnal experiment. 

(e) The relative dielectric constant of frozen ground kg Is 
assumed equal to 3. 

(f) The power transmission coefficient t_. at the snow-air 
boundary is calculated for each of the 21 measurements considered 
using. 




Cos 9-1 

Jk - Sln^ 9 


V s 


V^s 


Cos e + \/k^ + Sln^ e 


(4-15) 


where 


e = The look angle of the radiometer, 

= k' - jk'* = The comolex dielectric 
s s 


o'*. 20° and 50°. 
constant of snow. 
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is a function of the radiometer's look angle and the snow dielectric 
constant. The snow dielectric constant is calculated using the Tinga 
et al. (1973) mixing formula shown in Table 2-3, for each of the 21 
different measurements considered. is calculated as shown above 

So 

for the different values and are shown as a function of the 
percentage wetness by volume of snow in Figures 4-5, 4-6 and 4-7 for 
:? = 0°, 20° and 50° respectively. 

(g) The propagation angles in the snow layer are calculated us'ng 
Snell’s law given below as: 

ynr 

Sec e; = - 7 I— 5 — (4-16) 

''k^ - Sin"^ e 


The propagation an;le, 8 ^, is a function of the snow dielectric constant 
and, therefore, is a function of the snow wetness. 

(h) The power transmission coefficient at the ground-snow 
boundary is calculated using a similar expression to that shown by 
Equation (4-15), that is: 


gs 


1.0 - 


Cos 9 ; - \l 


'^^ks - Sin^ 9; 


^°^ ^s - Sin^ a; 


(4-17) 


kg is the dielectric constant of ground, x^^ is calculated for each 
of the 21 different measurement, jonsidered and for the different 
propagation angles in the snow layer. 

(i) The value of < is taken equal to that previously shown by 

d 

equation (4-13). 
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(j) The scattering coefficient is assumed to have the form: 


= A * 


(4-18) 


That is is assumed to have a linear dependence on the snow wetiess 
similar to < . A and B are constants to be determined. 

d 

(k) T^^, the downward emitted sky radiation scattered by the 

snowpack in the direction of the radiometer and previously shown by 
Equation (4-11), is assumed to be small in value and is neglected in 
this analysis. That is Tg, the brightness temperature of the snow and 
the ground underneath as described by Equation (4-10) is taken equal 
to/ T , the apoarent temoerature measured by the radiometer at e = 0*^, 
29° and 50°. This is because as shown in Figures 4-5, 4-6 and 4-7 
for as a function of at 3 = 0°, 20° and 50° respectively, the 
minimum value of x calculated is 9.9715 for 9 = 50° and M ('^) = 

S3 V 

2.53. This value when substituted into Eat :ion (4-11) and T^|^^ 
is taken equal to 300°K, the maximum possible value for T^^ 

is only equal to 8.55°K. This value is small enough to be neglected. 

T values measured are shown in Figure 3-1. 

3 p 

The values calculated by steps (a) - (k) above are substituted 

into Equation (4-14). This produces 21 equations, each describing the 

emission from the homogeneous layer of snow at a particular time through- 
<v 


out the diurnal experiment. These equations are solved using non- 
linear regression analysis ( BMDP , UCLA: 1977) to estimate the values 
of the constants A and B. This is done by solving Equation (4-14) 
for the eight measurements in which snow was totally dry in all its 
layers. In this case Equations (4-13) and (4-18) reduce to: 
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Figure 4-5 153 , the Power Transmission Coefficient at the Snow- 

Air Boundary as a Function of the Percentaqe 
Wetness by Volume of Snow at 6 = 0°. 
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0. 0.7 1.4 2.1 2.8 3.5 4.2 4.9 5.6 

Mv 1S> 

Figure 4-6 Tja* the Power Transmission Coefficient at the Snow- 
Air Boundary as a Function of My, the Percentage 
Wetness by Volume of Snow at 9 = 2QO. 
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Mv (•» 

Figure ^ • r^a* the Power Transmission Coefficient at the Snow- 
Air Boundary as a Function of My, the Percentage 
Wetness by Volume of Snow at 9 = 500. 
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and 


= 0.0203 

a 


Nep/cm 


(4-19) 


IC 


S 


= A 


, Nep/cm 


(4-20a) 


Once the value of A has been estimated, it is substituted into Equation 


(4-18) and the expression for with only constant B unknown is sub- 

£ 

stituted into Equation (4-14). Non-linear regression analyses are 

carried out for both the dry and wet snow conditions to estimate the 

e 

value of the Constant B. The results of these analy^s show that: 


A = 0.0171 Nep/cm (4-20b) 

B = 0 


The above equation gives a value for the scattering coefficient of 
snow at 37 GHz. tc^ predicted shows no dependence on snow wetness. 

Figure ^-8 shows a plot of k and < as predicted by this model 
as a function of the percentage wetness by volume of the snow layer. 

The values of k, and k. shown by Equations (4-13) and (4-20b) respec- 
tively are now substituted into Equation (4-10) and this equation is 
solved for the 21 different measurements to produce estimates of Tg 
at different times throughout the aiurnal experiment. Figures 4-9, 

4-10 and 4-11 show the observed T^p measurement along with the calculated 

Tg values as a function of time along the diurnal for 9 = 0®, 20® 

and 50® respectively. For all the three angles considered, the predicted 
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Figure 4-8 i<ra» <s ^nd as a Function of Snow Wetness as Predi<"* 
ted by the Single Homogeneous Snowlayer Model. 


76 


B 

§ 3.0 
i 2.0 
~> LO 
S 0 



Hour of Day 


Figure 4-9 Observed Tap and Tg Predicted by the Single Homogeneous 
Snowlayer Model as a Function of Time at 9 = 0°. 
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Figure 4-10 Observed T^p and Tg Predicted by the Single Homogeneous 
Snowlayer Model as a Function of Time at e - 200. 
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Figure 4-1’ Observed Tap and Tg Predicted by the Single Horuogeneous 
Snowlayer Model as a Functiorv of Time at e = 50°. 
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Tn response leads that of the observed T response when the snow is 
melting and lags behind it when the snow is refreezing. This is oje 
to the fact that the response of the predicted Tg is a function of the 
snow wetness of the top 5 cm layer of snow. The wetness of the top 5 cm 
layer of snov. is assumed here to be descriptive of the whole snowpack. 
However, *>e wetness of the too 5 cm layer <^f snow experimentally shows 
higher values of wetness than the rest of the sub-snow layers when the 
snow is melting and in turn shows lower wetness values than the rest 
of the sub-snow layers when the snow is refreezing. 

4 . 4 Emission ^tode1 for a Multilayered Snowpack 

As was previously mentioned in the last section, the snowpack 
during the 2/17 - 2/18/77 diurnal experiment had a total depth of 30 
cm and modeled to consist of four distinct snow layers. Figure 4-12 
shows the thickness of each of the snow layers within the snowpack. 

In this analysis. Equation (4-10) describing the emission from a 
multilayered snowpack is used to estimate the value of and its 
dependence on snow wetness. The snowpack in this case is treated 
as consisting of four distinct snow layers. 

Equation (4-10) is solved in a similar manner to that previously 
discussed in Section 4-3. In this case, however, the power transmission 
coefficients are calculated at every boundary between the snow layers, 
the boundary oetween the uppermost snow layer and the air, and the boun- 
dary between the ground and the lowermost snow layer. Similarly the 
propagation angles are calculated i>- each snow layer for the three dif- 
ferent look angles. Values of the thermometric temperatures in each 
layer and the wetness of each layer for the 21 different times 


80 




I 

I 

I 


Tc 

Radiometer 



Figure 4-12 A Geometrical Presentation of the Snowpack Cross Sec- 
tion During the 2/17/77 - 2/1C/77 Diurnal Experiment. 
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considered throughout the diurnal are shown in Table 5'1 md "inure 3-2 
respectively, o, the density of the snow was reoorted equal to 0.21 
q/cm^ for all the snow layers throuqhout the diurnal exoerinent. ’he 
results of this analysis gives 


: = 0.0175 - 0.0046 M 

S V 


'ieo/cn 


(4-21) 


That is the scatterinn coefficient, is smaller in value than the 
absorption coefficient, for drv snow. Also, increases linearly 

a o 

with increasing wetness but, decreases linearly with increasing wet- 
ness as shown in Figure 4 - 13 . The rate of decrease of with wetness 

is very small as shown by Equation (4-21). 

Figures 4-14, 4-15 and 4-16 show the T observed during the diurnal 

ap 

experiment alonn with the T^ predicted by the multilayer model at 

} =0°. and 50^ respectively. The values of Tg are calculated by 

substitutino tne values of and • of every layer into Equation (4-10) 

alonq with the other parameters to give the calculated values of Tg 

throuqhout the diurnal exoerinent. It would be seen from Figures 4-14, 

4-15 and 4-',5 tnat Tp calculated still leads T observed by the radiometer 

when the snow is -eltina and laas oehind it when the snow is refreezinq. 

However, the dif*er*ence in values between T observed and Tn predicted 

dp B 

is s'^aller than those previously shown for the single homogeneous layer 

model. The fits shown bv Figures d-i4, 4-15 and 4-16 are, therefore, 

better fits ^cr tl.ev produce estimates of Tg that are closer in value 

to the T "easuren hv the radiometer, 
ap 
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Figure 4-14 Observed Tap and Tg Predicted by the Multilayered 
Emission Model as a Function of Tine at 9 = 0°. 
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igure 4-16 Observed T^o and Tb Predicted by the Multilayered 
Emission Model as a Function of Time at 9 = 50°. 



4.5 Multilayered Model for Snow with Split Top Layer 

In order to modify the models applied in the previous two sections 
to give better estimates of the observed T^p, an understanding of the 
snow melting and refreezing process is necessary. Figures 3-3 to 3-6 
along with Table 3-1 indicate that as the air temperature above the 
snowpack increases sufficiently above freezing combined with the radia- 
tion from the sun shining on the snowpack, the snow begins to melt. 

The process of melting starts on the top of the snowpack with the top 
1 cm of snow or less becoming much more wet and so having a higher M^ 
than the rest of the snowpack below it. If the air temperature remains 
warm enough, the snow wetness increases with time in the top few 
centimeters of snow, becoming wettest on the top of the snowpack and de- 
creasing in wetness with increasing snow depth. At this time the radio- 
meter observing the snowpack at 37 GHz is sensitive to a small snow 
deptn on the top of the snowpack. Due to the decrease in penetration 
depth, this very small snow depth is responsible for the majority of 
the emission of microwaves from the snowpack. This top and very 
wet snow layer could be smaller in thickness than the top 5 cm layer 
of snow considered in the previous model and whose average wetness 
was measured by Stiles and Ulaby (1980) and is shown in Figure 3-2. 

As the air temperature starts decreasing later in the afternoon 
to a temperature below that of freezing, the snowpack starts to refreeze 
The refreezing process, similar to the melting process, starts from the 
top of the snowpack. A small depth of snow on the top of the snowpack, 
which could be 1 cm in depth or smaller, starts to refreeze and the 


87 



of that depth decreases greatly. With time, the wetness, M^, of this 
small depth drops to values that are less than the values of of the 
lower depths of the snowpack. The snow in the top of the snowpack then 
becomes dry and its drops to zero in the top few centimeters of the 
snowoack, while the snow in the lower layers of the snowpack remains 
higher in wetness. With time and continuing decrease in air temperature, 
the frozen layer on the top of the snowpack increases in thickness 
gradually down the snowpack until the whole snowpack is frozen. 

The radiometer observing the snowpack at the start of the refreezing 

process will be sensitive at times to a snow depth that is larger than 

the 5 cm top layer of snow due to the increased penetration depth of the 

wave with the disaopearance of free water from the top of the snowpack. 

The snow wetness measurements collected from the top 5 cm layer of snow 

are average values of wetness throughout the snow thickness considered 

and do not give the exact wetness profile in the snow thickness. This 

average value of wetness will, therefore, lag behind the T^p observed 

when the snow ’s molting and lead the T observed when the snow is 

ap 

refreezing as shown in Figure 3-1. This causes the hysteresis pattern 
shown in Figures 3-7, 3-8 and 3-9. 

The node! Jevelooed as described by Equation (4-10) needs to be modi- 
fied to account for the ohysical behavior of snow at its nrelting and 
refreezinc oracess. This is done hy splitting the uopermost snow layer, 
shown in Fiaure 4-1^ as Layer 4, into two sublayers each 2.5 cm in 
thickness. The thermometric temperature of these two sublayers is 
taken to be the same as the origina 5 cm top layer. The percentage 
wetness by volume of the top 2.5 cm layer, shown in Figure 4-17 
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Figure 4-17 A Geometrical Presentation of the Revised Snowpack Cross 
Section During the 2/17/77 - 2/18/77 Diurnal Experiment. 

The Top 5 cm Layer is Shown Split into Two 2.5 cm Sublayers. 
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as Layer 5, is made seven times as wet as the original 5 cm top layer 
in the melting stage and then allowed to saturate at a value of = 

2.53 which was the maximum value of reached by the original 5 cm top 
layer. The percentage wetness of the lower 2.5 cm layer, shown in 
Figure 4-17 as Layer 4, is calculated during the melting process such 
that the average value of wetness of both 2.5 cm snow layers is 
equal to that of the original 5 cm uppermost snow layer. The wetness 
of the second 2.5 cm layer is made seven times as wet as the original 
5 cm uppermost snow layer during the refreezinq stage allowing it to 
saturate at the same value of = 2.53. The wetness of the top J.5 
cri layer during the refreezing stage is then calculated such that 
the average wetness of both 2.5 cm snow layers is equal to that of 
the original 5 cm uppermost snow layer. Figure 4-17 shows the snow layers 
within the snowpack and their various thicknesses. Figure 4-18 shows 
the calculated wetness curves of the two 2.5 cm layers that make up the 
original 5 cm uppermost snow layer, along with the wetness curve of the 
original 5 cm top layer of snow. 

Equation (4-10) describing the emission f : im a multilayered snowpack 
is solved in a manner similar to that shown in the last two sections, 
taking into account the modifications introduced in this section. The 
scattering coefficient estimated by this model as modified above is the 
same as that estimated in the previous section and shown by Equation (4-21). 
Figures '19, 4-20 and 4-21 show the T^p observed during the diurnal 
experiment along with the Tg predicted by the multilayered model with the 
split top layer at 9 = 0®, 20° and 50° respectively. The difference 
between the ^its shown here by Figures 4-19, 4-20 and 4-21 and those 
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Figure 4-18 Percentage Wetness by Volume My Curves for the Top 
2.5 cm Snow Layers Along With the Original 5 cm 
Snow Layer Wetness Curve of the 2/17/77 - 2/18/77 
Diurnal Experiment. The Top 2.6 cm Layer is Seven 
Times as Wet as the Original 5 cm Layer in the Melting 
Stage While the Lower 2.5 cm Layer is Seven Times as 
Wet as tiie Original 5 cm Layer in the Refreezing Stage. 
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Figure 4-19 The Observed Tap and the Predicted Tg as Calculated 
by the Multilayered Emission Model With Split Top 
Layer at 0 = OO. The Uppermost Sublayer is 2.5 cm 
Thick and Seven Times as Wet as the Original 5 cm 
Top Snow Layer. 
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Figure 4-20 The Observed Tap and the Predicted Tb as Calculated 
by the Multilayered Emission Model with Split Top 
Layer at e = 20°. The Uppermost Sublayer is 2.5 cm 
Thick and Seven Times as Wet as the Original 5 cm 
Top Snowlayer 
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Figure 4-21 The Observed Tap and the Predicted Tb as Calculated 
by the Multilayered Emission Model with Split Top 
Layer at 9 = 50O. The Uppermost Sublayer is 2.5 cm 
Thick and Seven Times as Wet as the Original 5 cm TOj. 
Snowlayer. 
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calculated previously and shown by Figures 4-14, 4-15 and 4-16, for 
e = 0*^, 20® and 50® respectively, for Tg are very small. 

The uppermost 2.5 cm snow layer was made eight, 12 and 15 times 
as wet as the original 5 cm uppermost layer of snow. The original 
top 5 cm layer of snow was split into two sublayers. The thickness 
of the top layer was made 1 cm and that of the second layer 4 cm in 
thickness and the analysis repeated again for the different wetnesses 
considered. However, the fits produced by the different thickness 
of the snow sub-layers or the different wetnesses considered were not 
much different from those shown by Figures 4-19, 4-20 and 4-21 and are, 
therefore, not shown here. 
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5.0 EVALUATION OF THE EMISSION MODELS AND THEIR RESULTS 


In this chapter, all of the emission models developed and applied 
in Chapter 4 to describe the emission from a snowpack are evaluated. 

The model that predicts the best estimates of the T observed by the 

dp 

37 GHz radiometer throughout the diurnal experiment is considered the 
superior model. The values of the absorption coefficient and 

a 

ic^, the scattering coefficient estimated by the superior model chosen 
are comoared to similar values reported in the literature and are 
used to calculate the percentage contribution from each layer of the 
snowpack throughout the diurnal experiment. 

5.1 Evaluation of the Models 

Two methods are used to evaluate the models developed in Chapter 

4. The first is a residual sum of squares analysis in which the sum 

of the squares of residuals occuring between Tg(N+) predicted by the 

model and T observed is calculated for each model. The model that 
ap 

shows the least residual sum cf squares at all of the radiometer's 

look angles is therefore considered a better model for it gives the 

best estimates of T,„ observed. Table 5-1 shows the results of such 

ap 

an analysis conducted on all the models developed. The model describ- 
ing a multilayered snowpack with split uppermost layer of snow and 
whose 2.5 cm top layer is seven times as wet as the original 5 cm 
uppermost layer of snow shows the least residual sum of squares for 
all 

The second evaluation method conducted, is a linear correlation 
analysis between the values of the Tg(.N^) predicted by each model and 
the T observed by the radiometer. Figures 5-1 to 5-9 show the results 

d D 
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TABLE 5-1 


The Residual Sum of Squares Analysis 
Conducted for all the Models Developed and 
at 0 = QO, 200 and 5QO. 


Model 

OO 

200 

o 

o 

in 

Single homogenous layer 
Snowpack 

2044.23 

3732.51 

9022.11 

Multilayered Snowpack 

1690.29 

2077.31 

3220.84 

Multilayered Snowpack 
with Split Top Layer 

1380.06 

1598.42 

3148.98 1 

i 

■ ■ > j 
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T, (Prad) ■ 1.010 i;,p(Obs) + 14.751 



Observed Tap in °K 


Figure 5-3 Linear Correlation Coefficient Between the Predicted 

Tb as Calculated by the Single Honwgeneous Layer 

Snov^pack Model at e = 50° and the T Observed. 

ap 
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To (Prso) ' 0.979 Top(Obs) + 3,692 



150 170 190 210 230 250 270 290 
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Figure 5-4 Linear Correlation Coefficient Between the Predicted 

Tb as Calculated by the Multilayered Snowpack Model 

at e = 0° and the Observed. 

ap 
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Tb (Prso)- 0.979 Tap(Obs) + 3, 004 




170 


_j I I - I . 

190 210 230 250 

Observed Tap in ®K 


-1 ! 

270 290 


Figure 5-5 Linear Correlation Coefficient Between the Predicted 

Tb as Calculated by the Multilayered Snowpack Model 

at 0 = 200 and the Observed. 

ap 
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I’e (Praci) • 0.987 T„p(Obs) + 11.640 



Figure 5-6 Linear Correlation Coefficient Between the Predicted 

Tb as Calculated by the Multilayered Snowpack Model 

at 9 = 50° and the T ^ Observed. 

ap 
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T, (ProG) • 0.970 Top(Obs) +5.560 



150 170 190 210 230 250 270 290 


Observed Tap in *K 


Figure 5-7 Linear Correlation Coefficient Between the Predicted 
Tb as Calculated by the Multilayered Snowpack with the 
Split Top Layer and the T^g Observed at e * 0°. New 
Uppermost Layer is 2.5 cm in Thickness and Seven Times 
as Wet as the Original 5 cm Top Snow Layer. 
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Tj, (Prso) • 0.957 Top(Obs) +9.215 



150 170 190 210 230 250 270 290 

Observed Topin ®K 


Figure 5-8 Linear Correlation Coefficient Between the Predicted 
Tg as Calculated by the Multilayered Snowpack with the 
Split Top Layer and the Tgp Observed at 0 = 20 ^. 

New Uppermost Layer is 2.5 cm in Thickness and Seven 
Times as Wet as the Original 5 cm Top Snow Layer. 
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Tb (Prao) « 0.966 Top(Obs) + 16.494 



Figure 5-9 


Linear Correlation Coefficient Between the Predicted 
T 3 as Calculated by the Multilayered Snowpack with the 
Split Top Layer and the T^p Observed at e = 50°. 

New Uppermost Layer is 2.y cm in Thickness and Seven 
Times as Wet as the Original 5 cm Top Snow Layer. 
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of such an analysis carried out on all the results of the models devel- 
oped at 0 = 0°, 20° and 50°. Figures 5-1, 5-2 and 5-3 show the linear 
correlation between the Tg predicted by the single homogeneous layered 

snowpack and the T „ observed at 0 = 0°, 20° and 50° respectively, 
ap 

Figures 5-4, 5-5 and 5-6 show the linear correlation beU-een the Tg 

predicted by the multilayered snowpack and the T observed at e = 0°, 

ap 

20° and 50° respectively. Figures 5-7, 5-8 and 5-9 show the linear 

correlation between the Tg, predicted by the multilayered snowpack 

model with the original uppermost snow layer split into two sublayers, 

each 2.5 cm in thickness, and the T,„ observed at 0 = 0°, 20° and 

ap 

50° respectively. Both multilayered emission models show superior 
estimates of Tg than the single homogeneous layer model. However, there 
is little difference in the correlation coefficient between the values 
of To predicted by both multilayered models and the T^^ observed. 

u ap 

On this basis the multilayered snowpack model, with the uppermost layer 
split into two sublayers each 2.5 cm in thickness and the top sublayer 
is made seven times as wet as the original 5 cm uppermost snowlayer 
is considered the best model of the three developed. 

5.2 Evaluation of the Values of k and Estimated by the Multilayered 
a s 

Snowpack Model 

Equations 4-13 and 4-21 show the absorption coefficient calculated 

by the Tinga et al. (1973) mixing formula and the scattering coefficient 

estimated by the multilayered emission model of snow as a function of 

snow wetness. The values of k,, and calculated in this study 

a s e ^ 

for dry snow are compared to those previously reported in the literature 
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as shown in Table 5-2. The values of k., k and calculated for dry snow 
in this study are greater in value than those reported by both Stiles 
and Ulaby (1980) and Hofer and Hatzler (1%0). This is due to the following 
reasons: 

(a) The difference in the snow deptii, snow density, snow t^|»rature, 
and snow crystal structure in all the experiments conducted. There 

were also differences in air temperature and other environmental con- 
ditions which could have caused the differences seen in the vaicns 
of and Kg. 

(b) The values of k., k_ and k^ estimated by Stiles and Ulaby 

o S C 

(1980) are obtained from the results of the snow pile experiment shwn 

in Figure 3-10. In their experiment, the snow was piled up to a depth 

of 170 cm and the apparent tonperature of tha snow pile at 37 GHz 

war measured for different depths. At the same time, T^^ was measured 

fr ...I the undisturbed snowpack around the snow pile to see the difference 

in the T.„ measurement for the same depth. Figure 3-10 shows that for 
«P 

a snow depth, above the ground, of 30 cm or for a 6.3 cm water equivalent 

at e = 27°, the T measured from the snow pile was 260°K whereas the 
ap 

T measured from the undisturbed snow was 200°K. Similarly for 
ap 

0 = 57° and a water equivalent of 15 cm, T^p measured from the snow 

pile was 200°K whereas T__ measured from the undi studied snow was 175°K. 

ap 

That is the values of T measured from the snow pile are higher than 

ap 

those measured from the undisturbed snow. This difference is due to 
the alteration in the crystal structure of the snow as it was piled 
causing the scattering coefficien .o become smaller in the snow pile 
than it is for the undisturbed snow and, t*- * ^'v;ng higher 
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TABLE 5-2 


Values of K'a, •:$. <e for Dry Snow 

as a Function of e Calculated Usinq Stiles and Ulaby (1980) Model, 
the Multilayered Snowpack Emission Model 
and the Hofer and Matzler (1980) Model. 


Model 

4 

In Degrees 

k* 

a 

dB/cm 

■'s 

dB/cm 

•'a'^'e 

*^e 

dB/cm 

Stiles and Ulaby (1980) 

27® 

0.04 

0.04 

0.52 

0.08 

Stiles and Ulaby (1980) 

57® 

0.07 

0.05 

0.59 

0.12 

Multilayered Emission Model 

0®, 20® & 50® 

0.18 

0.15 

0.55 

0.33 

Hofer & Matzler (1980) 
(High Winter) 

0®-60® 

0.02 

0.03 i 

1 

0.40 

0.05 

Hofer & Matzler (1980) 
(Spring) 

0®-60® i 

( 

0.02 

0.04 1 

0.33 

0.06 

1 
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values of T.^ from the snow pile than undisturbed snow. Since the 
ap 

values of k , k and k calculated in this study are those of undis- 

G S 6 

turbed snow while those calculated by Stiles and Ulaby (1980) are those 
for the snow pile* the larger values of and < calculated here 

are, therefore, explainable. Figures 5-10 and 5-11 show the brightness 
temperature calculated using Equation (4-14) for a snowpack made up 
of one homogeneous layer of snow as a function of water equivalent 
at 0 * 27° and 57° respectively. These figures are calculated by 
substituting the values of and obtained in this study for 

undisturbed snow into Equation (4-14). The figures also show the measure- 
ments of Stiles and Ulaby (1980) for undisturbed snow at 9 = 27° and 
57° of Figure 3-10. Notice the snail difference in Tg calculated here 
and that measured by Stiles and llaby (1980) at the two angles con- 
sidered. 

(c) The ratio of k ./ k . is very similar in all the measurements 

a c 

shown in Table 5-2, however the value of tends to be much different. 

The dependence of Tg on is shown in Figure 5-12. In this figure, 

Tg is plotted versus for a snowpack made up of one homogeneous 

layer and 30 cm in depth and a snow temperature of 262.33°K. k 

calculated in the study is comparable in value to the total loss measured 

by Stiles and Ulaby (1980) for dry snow as shown in Figure 2-10. This 

adds to the credibility of the value of calculated here. 

e 

5.3 Total Emission from Each Snow Layer 

The emission from each snow layer within the snowpack, shown 
in Figure 4-9 is calculated using the values of and k. given by 

d S 


no 




Water Equivalent, W in cm 


Figure 5-10 To Calculated From an Undisturbed Snowpack as a Function 
OT W for One Layered Snowpack at 6 = 27° Using the 
<a. <s afid Calculated in This Study. 
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Figure 5-11 To Calculated From an Undisturbed Snowpack as a Function 
of W for a Single Layer Snowpack at 0 = 57® Using the 
<3. Kj and Calculated in the Study. 


112 




Figure 5-12 Tg as a Function of <e ® Snowpack of One Homogen- 
eous Layer at Nadir. The Snow Depth is 30 cm and Snow 
Temperature is 262.330K. 
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Equations 4-13 and 4-21 respectively, and Equation (4-10). Figure 
5-13 shows the percentage contribution of each layer within the snow- 
pack as a function of time throughout the di»rnal experiment. It is 
seen that for dry snow, the ground contributes the most to the total 
microwave emission measured by the radiometer. The contribution from 
the ground underneath a totally dry snowpack is above 45% of the total 
emission measured. As the topmost snow layer of the snowpack increases 
in wetness, its percentage contribution to the total emission increases 
as well. When the wetness of the top 5 cm layer is 1% by volume it 
contributes almost 80% of the total emission. When the snow wetn ss 
of the too 5 cm layer of snow is 2% by volume, it contributes 90% 
of the total emission. For dry snow, therefore, the emission from the 
snowpack and the ground underneath it is a volume emission stemming 
from each layer within the snowpack and the ground underneath it. 
However, this volume emission changes rapidly with the appearance of 
free water in the uppermost layer of the snowpack to become a surface 
emission stemming almost totally from the surface layer of the snowpack. 
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Figure 5-13 The Percentage Contribution of Each Layer Within the 
SnoMpack as a Function of Time. 
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6.0 CONCLUDING REMARKS 


In this paper a model describing the microwave emission from 
snow at 37 GHz was developed as a function of the physical and dielectric 
parameters of the snowpack. These parameters included the consideration 
of a snowpack consisting of different layers. The model provided 
information on the contribution of each snow layer and the ground to 
the total emission, for both wet and dry snow conditions. The model 
also provided solutions of the scattering, absorption and extinction 
coefficients for both wet and dry snow. The errors that exist in this 
model are due to the following: 

(a) The nndel develooed predicted the brightness temperature 
of the snow and not the apparent temperature actually measured. 

It, therefore, did not take into account the contribution of ttie downward 
emitted sky radiation scattered by the scene in the direction of the 
radiometer antenna. 

(b) The lack of information of the dielectric constant of wet 
snow at 37 GHz, and of the dielectric constant of each snow layer 
considered throughout the experiment. 

(c) The lack of wetness profiles that provide the wetness of 
the snowpack accurately and for every 1 cm thickness of the snowpack. 

(d) Lack of information on the exact contribution of the ground 
to the apparent temperature measured by radiometer, since no T 

ap 

measurements were conducted for the frozen ground of the scene during 
the experiment. 
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The model is the first to calculate values of < , < and k as a function 

a S 0 

of snow wetness. The model is also the first to show the contribution 
of each snow layer to the total emission from snow and ground as a 
function of snow wetness. 
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